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We use angle-resolved photoemission spectroscopy (ARPES) to image the emergence of substaintial disper- 
sion anomalies in the electronic renormalization of the actinide compound UCoGas which was presumed to 
belong to a conventional Fermi liquid family. Kinks or abrupt breaks in the slope of the quasiparticle dispersion 
are detected both at low (~130 meV) and high (~1 eV) binding energies below the Fermi energy, ruling out 
any significant contribution of phonons. We perform numerical calculations to demonstrate that the anoma- 
lies are adequately described by coupling between itinerant fermions and spin fluctuations arising from the 
particle-hole continuum of the spin-orbit split 5/ states of uranium. These anomalies are resemble the 'wa- 
terfall' phenomenon of the high-temperature copper-oxide superconductors, suggesting that spin fluctuations 
are a generic route toward multiform electronic phases in correlated materials as different as high-temperature 
superconductors and actinides. 



I. INTRODUCTION 

The coupling between electrons and elementary excitations, 
originating from lattice or electronic degrees of freedom, can 
drive the formation of new emergent phases, such as mag- 
netism and superconductivity. For rare-earth and actinide f- 
electron systems, the development of low-energy fermionic 
excitations with heavy electron mass is complicated by the 
interactions between / electron spins and those of itinerant 
electrons. They often exhibit an interplay between Kondo 
physics, magnetism, and superconductivity Q. So far it has 
been thought that this physics is not at play in UCoGa 5 with 
a tetragonal crystal structure and metallic ground state Q. In 
fact, it has been dubbed a 'vegetable', because of its lack of 
heavy mass (3] 01, other competing orders 0EHS1, and show- 
ing of a spin-lattice relaxation rate that obeys the Korringa law 
of a conventional Fermi liquid at low temperatures |8|. How- 
ever, about ~ 75 K it violates the Korringa law by exhibit- 
ing a quadratic temperature dependence without any sign of 
magnetic ordering. The simple lack of ordering is often con- 
trasted with the rich properties of the isostructural PuMGas 
and PuCoIn5 compounds (M=Co, Rh) [9 ]. Here, we show that 
UCoGas is not a typical Fermi liquid. It is rather anomalous 
as evidenced by the drastic renormalization of the electronic 
dispersion at low and high binding energies when compared 
with ab-initio electronic bands. These kinks resemble the 'wa- 
terfall' phenomenon of cuprates [11, d[l8) and point toward 
a common mechanism of spin fluctuations in both d and / 
electron systems. 

The actinide systems remain poised between the strong and 
weak coupling limit of Coulomb interaction. Therefore they 
offer significant tunability across several correlated states of 
matter. The effective Coulomb repulsion U of the 5/ elec- 
trons is not strong enough to localize all of them, yet, it is 
sufficient to slow them down so they acquire a moderately in- 
creased mass near the Fermi level. The central puzzle is how 
does this common low-lying heavy electronic state transform 
into stable ground states which differ dramatically within the 
family (6l[T3]|. An understanding of these differences hinges 
on knowing the nature of the exchanged bosons that dress the 
bare electrons to become heavy quasiparticles and ultimately 
drive the system into a particular stable ground state. The 



choice of UCoGas for this study allows us to rule out any pos- 
sible intervention of competing orders. Thus it offers a clean 
approach to unravel the nature of the exchanged boson in the 
actinide family and related heavy-fermion systems. 



II. EXPERIMENT 

We use ARPES to image the momentum and energy de- 
pendence of the electronic dispersions and thereby reveal the 
renormalization of itinerant bands. Details of the ARPES ex- 
periments and calculations are described in the supplementary 
material (SM |[T4ll ). Figure 1 reports our main experimental 
result of UCoGas . A quick visual inspection of the spectral 
function in Fig. [TJa) and [TJb), (3D and 2D rendering of the 
same data, respectively), shows extensive spectral weight re- 
distribution both in momentum and energy space. The dis- 
persion (traces of the intensity peaks), in Fig. []Jb), reveals a 
drastic departure of the quasiparticle states from the ab-initio 
electronic structure calculations (black and dashed). More im- 
portantly, the associated quasiparticle width at the peak posi- 
tions is significantly momentum and energy dependent. Un- 
like the renormalized quasiparticle dispersion relation, which 
is a genuine manifestation of many-body correlation effects, 
the anomaly in the lifetime may have many intrinsic and ex- 
trinsic origins and is more difficult to assign to a particular 
mechanism [ 1 ] [see SM[ 14] on the extrinsic background]. We 
notice that while bare bands, computed within generalized 
gradient approximation (GGA) based density functional the- 
ory, are present in the entire energy and momentum region, the 
spectral intensity is accumulated mainly at two energy scales. 
Such a connection between low and high energy scales is anal- 
ogous to the so-called 'waterfall' or high-energy-kink feature 
observed in single-band cuprates [Gj][l2l[T8]], where the band 
bottom lies at T, while here it is shifted to the M point. 

The anomaly is markedly different in the energy and mo- 
mentum space, which is a hallmark feature of correlated 
electron states. To expose this distinction, we present sev- 
eral energy-distribution curves (EDCs) of the ARPES inten- 
sity at several representative fixed momenta in Fig. []Jc), and 
momentum-distribution curves (MDCs) at several fixed en- 
ergy points in Fig.[TJd). The peaks in the EDCs display promi- 
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FIG. 1. Dispersion and measured ARPES spectral function anomaly, (a) 3D intensity map of UCoGas along the M— ^ F direction (h, h, 0) 
in the Brillouin zone, (b) Same data plotted as 2D contour map and compared to the corresponding ab-initio GGA electronic band- structure 
dispersions (black dashed lines). Arrows indictate two quasi-non-dispersive energy scales of higher intensities. All the ARPES data are 
collected at the photon energy of 46 eV, with energy resolution better than 20 meV. (c) The EDCs are plotted for fixed momenta. The curves 
from bottom to top are chosen for equally spaced momentum points from M to F. Arrows have same meaning as in (a), (d) MDCs for fixed 
binding energy. Bottom to top curves are chosen from E=-2.1 eV to E=0 eV. The dashed lines are guide to the eyes for two extracted low- 
lying dispersions of anomalous character, analyzed in Fig. 2. (e) Photoemission spectrum (PES) of UCoGas at photon energies his=\08 eV 
(red) and hv-Wl eV (blue). Their difference curve (green) reveals the dominant 5/ character 1 1 1. (f) Partial PES of UCoGas (momentum 
integrated between F and M) compared with two related actinides. UNiGas and PuCoGas exhibit antiferromagnetic [13 1 and superconducting 
1 18 1 ground states, respectively. All data are normalized arbitrarily for ease of comparison. 



nent dispersionless features at two energy scales. The lowest 
energy peak persists at all momenta around -80 meV and grad- 
ually becomes sharper close to the F point. This low-energy 
feature reveals the formation of the long-lived renormalized 
quasiparticle. We will show below that the renormalization 
phenomena can be quantified to originate from the interac- 
tion of electrons with spin fluctuations, even without invoking 
the more traditional Kondo physics of heavy fermions. The 
high-energy feature around —1 eV is considerably broad in 
both energy and momentum space. In momentum space it 
attains a large width, demonstrating that these states are sig- 
nificantly incoherent. By comparing the experimental dis- 
persion with its ab-initio counterpart, we can convincingly 
draw the conclusion that these incoherent states are created 
by spectral- weight depletion near —500 meV due to coupling 
to spin-fluctuations. To gain further confidence about the cor- 
relation origin of these band anomalies, we note that such a 
redistribution of the spectral weight does not occur uniformly 



for the non-interacting bands; for instance the spectral weight 
is strongly suppressed between —2 eV to —1.5 eV near the 
M point, albeit several bands are present in this part of the 
Brillouin zone. By comparing the MDCs and EDCs, we im- 
mediately see that the spectrum is less dispersive as a function 
of momentum, while it disperses strongly with energy. 

Figure [TJe) shows the quasiparticle peak pushed below the 
Fermi level by approximately 290 meV in the integrated pho- 
toemission spectrum (PES). The difference curve (green line) 
between measurements at 108 eV and 102 eV photon energies 
reveals that the quasiparticle states in the vicinity of the Fermi 
level are predominantly of 5/ character (2. The partial PES 
that comes from the F — »M momentum direction, is obtained 
by integrating the spectrum in Fig.[TJb) and shown by the red 
line in Fig.[TJf). Interestingly, it displays a characteristic peak- 
dip-hump feature, which is generic for other isostructural ac- 
tinides (6l[T6) (also shown here), as well as for £-Pu ifTTl and 
copper-oxide superconductors (TT1 [T2j |T8]| (not shown). Un- 
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FIG. 2. The real and imaginary parts of the experimental self-energy, 
(a) Two low-lying dispersions derived from the peak positions of the 
MDCs [Fig.fTJc)] along M— ^ T. The momentum axis is offset by 
for band Bl. Dashed lines of same color give the ab-initio disper- 
sion (offset from their corresponding Fermi momenta values) which 
help expose the degree of anomaly for each dispersion with respect 
to their linear dispersion features, (b) The real part of the self-energy 
measures the deviation of each dispersion curve in (a) from its cor- 
responding GGA value. Green arrow marks the location of the low- 
energy kink, (c) The inverse of the MDCs width obtained with the 
help of Fig. [T| C X which is proportional to the imaginary part of the 
self-energy, when the extrinsic background contribution is negligible. 
Different color fillings in all three curves (corresponding to band Bl) 
separate the characteristic energy scale (marked by red arrows). The 
imaginary part attains a peak at the energy where the real part of the 
self-energy changes sign. To demonstrate this behavior we calculate 
Im E from Re E in (b) and plot is as black dashed line in (c). 



like UCoGas, the actinides UNiGas and PuCoGas exhibit an- 
tiferromagnetic fT3l and superconducting [18 ] ground states, 
respectively. 



III. ANALYSIS 

In order to explain the correlation aspects of the dispersion 
anomalies, we extract two quasiparticle dispersion branches 
by tracing the locii of the MDC peaks as shown in Fig. [2ja). 
To disentangle subtle features in the dispersion, the common 
practice is to observe the departure of the dispersion with re- 
spect to a featureless one, £ n fc (shown by dashed lines of same 
color). A visual comparison reveals a sharp kink in the metal- 
lic actinide UCoGa 5 , that is, a change in the slope of the dis- 
persion. The low-energy kink lies around —130 meV (indi- 
cated by green arrow), which is higher than that in cuprates 
(—70 meV) H9l . On the other hand, the kink observed in the 
actinide USb2 is at much lower energy, and thus is most likely 
caused by electron-phonon coupling [20]. Furthermore, the 
high-energy kink of the 'waterfall' (—1 eV) poses a ubiquitous 
anomaly as in cuprates (TTl[T2j[T8]|. On the basis of these com- 
parisons, we deduce that the present kink lies well above the 
phonon energy scale of ~30 meV for metallic UCoGas EL 



and hitherto provides a novel platform to study the evolution 
of spin fluctuations in /-electron systems. It it known that the 
role of phonon vanishes in the high-energy scales of present 
interest in conventional metals, unlike, for example, in quasi- 
one-dimensional blue and red bronzes where Peierls instabil- 
ity can shift the phonon energy scale by 10-15 times. l2Tll22ll 
In standard quantum field theory notation, the coupling be- 
tween bosonic excitations and fermionic quasiparticles is de- 
fined by a complex fermionic self-energy E n (fc,o;) function. 
The real and imaginary parts of the self-energy cause dis- 
persion renormalization and quasiparticle lifetime broaden- 
ing, respectively. This information can be derived from the 
ARPES spectra by tracing the peak positions of the n-th band 
kfi as 

Re E n (fc, lj) = uj - £ k ,n (1) 
ImE n (fe,w) oc 1/FWHMmdc- (2) 

Here £fc jn is the GGA dispersion with respect to its corre- 
sponding Fermi energy. The obtained results for Re E and 
the inverse full width at half maximum (FWHM) of MDC are 
given in Fig.[2jb) and[2tc), respectively, for the two extracted 
bands shown in Fig. Qa). Taking advantage of our high- 
precision data, we expose two energy scales in Re E (marked 
by arrows) which are intrinsically linked to the quasiparticle 
lifetime r = k/ (21m E). In addition to the peak at the kink 
energy, Re E possesses a sign reversal around 260 meV. This 
is an important feature which imposes the constraint that the 
corresponding Im E should yield a peak exactly at the same 
energy, due to the Kramers-Kronig relationship. This is in- 
deed in agreement with our findings for UCoGas in Fig.[2jc), 
which exhibits a sharp peak in the extracted broadening, ex- 
actly where Re E changes sign (marked by green arrow). It is 
worthwhile to mention that in ARPES measurements, the ex- 
trinsic source of quasiparticle broadening is typically known 
to be quasi-linear in energy [24], while the presence of a peak 
is a definitive signature of the intrinsic origin of broadening. 



IV. NUMERICAL CALCULATIONS 

Encouraged by the aforementioned analysis, we proceed 
with an ab-initio calculation to delineate the nature of the el- 
ementary excitations responsible for this anomaly. As men- 
tioned before, both the energy scales and the strength of 
electron-phonon interaction (121 Q2) is insufficient to capture 
our observed dispersion anomalies. We demonstrate that the 
electron-electron interaction due to dynamical spin fluctua- 
tions of exchanged bosons can describe the data. In UCoGa 5 
strong spin-orbit coupling from relativistic effects enables 
substaintial spin fluctuations whose feedback results in an in- 
crease of the electron mass and a shortening of the quasipar- 
ticle lifetime. We perform numerical calculations of the spin- 
fluctuation interaction potential (defined here as V e f /) and re- 
sulting self-energy E n (fe,o;) within spin-fluctuation coupling 
theory. Our intermediate coupling approach constrains the 
value of the Coulomb potential U to be of the order of the 
ab-inito band width ~ 1 eV, for details see lfT4l . 
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FIG. 3. Theoretical study of the spin-fluctuation dressed dispersion anomaly, (a) Spectral function intensity map versus energy and momentum 
along M— ^ T. Black lines are the corresponding ab-initio bands. Arrows mark the two energy scales of highest intensity, as in Fig. [l|a). 
(b) The real part of the momentum and orbital averaged self-energy (blue line), compared with the experimental data (red-dashed line). The 
experimental curve is multiplied by a factor of 5 for better visualization, because it is derived from a linear slope, (c) The imaginary part 
of the self-energy shows a peak at the same energy where the experimental linewidth also has a peak, (d) Computed partial PES (integrated 
over momentum between V and M) shows multiple peaks, as in the experimental curve (red dashed lines, taken from Fig.[T£). The black line 
is the corresponding (normalized) GGA density of states. The red arrows mark the depletion of states, originating from the spin fluctuation 
interaction. 



The computed results are presented in Figs.[3jb) and[3jc). 
A close examination of all these spectra reveals essentially 
two energy scales tied to the electronic structure of UCoGa 5 . 
In this system, the spin-orbit coupling of U atoms splits the de- 
generacy in the 5/ states by ~ 500 meV |9 ], which stipulates a 
strong peak in V e ff, and Im E (correspondingly Re E changes 
sign at this energy owing to the Kramers-Kronig relation). 
The second peak develops at a higher energy from the tran- 
sition between unfilled 5/ states and occupied 3d states of the 
Co atoms. Note that both energy scales are clearly discernible 
in V e f f (red arrows in the insert) and shift to higher energy 
in ImE due to the effects of dispersive bands. The peaks 
(dips) in Im E correspond to electronic spectral- weight dips 
(peaks) in ARPES as shown in Fig.[3ja). Importantly, the spin 
fluctuations also adequately explain the quasi-dispersionless 
renormalized feature around -80 meV and the spectral- weight 
peak-dip-hump features (marked by arrows), as observed in 
the ARPES spectra in Fig. [TJb), although some discrepan- 
cies exist. The momentum-integrated spectrum in Fig. [3Jd) 
further demonstrates the presence of both dips in the exper- 
imental and theoretical spectra. We note that the calculated 
self-energy corrections only shift the GGA Fermi surface (FS) 
maps in the k z = plane. The computed FS maps are in ac- 
cord with de Haas- van Alphen (dHvA) measurements, except 
for the tube-like FS around (7r/2, 7t/2, 0) lfT6l . Furthermore, 
our GGA calculation agrees with previous electronic structure 
calculations (ElE), except for the detailed shape of the tube- 
like FS. By turning on the self-energy the disconnected tube- 
like FSs become connected along the k z direction, as shown 



in the SM [14], signaling a topological transition between a 
closed and open FS due to a small shift in Fermi energy. 

Next we turn to estimating the degree of correlations in this 
material. The electron-phonon and spin- fluctuation coupling 
strength A can be extracted from the slope of the self-energy 
in the low-energy region as Re E n (fe, £n(fcp)) = — A£ n (fcp). 
Our theoretical estimate of the spin- fluctuation coupling con- 
stant is Xsf = 1-0. Experimentally, an estimate of A at 
the Fermi level is obtained by comparing the measured Som- 
merfeld coefficient with respect to its non-interacting value, 
7 « (1 + A)7o. The reported measured values of 7 are 
between 7 =10-21 mJ/mol K 2 0. Our GGA calcu- 
lated value of 7o=5.8 mJ/mol K 2 agrees well with the lit- 
erature (l2l Q31 and thus gives an experimental estimate of 
X exp ^0.7 — 2.6. Therefore, the combination of spin fluctua- 
tions, Xsf = 1-0, and low-energy electron-phonon coupling, 
^ph ^0-7 for PuCoGa 5 (23H28), adequately describes the ob- 
served electronic mass renormalization in UCoGas . 



V. CONCLUSIONS 

Our intermediate Coulomb-coupling approach offers an op- 
portunity to address the complex nature of the dynamic band 
renormalization in actinides. This is not unexpected for ma- 
terials with itinerant electrons, since similar theories based on 
the coupling between spin fluctuations and fermionic quasi- 
particles are among the leading contenders for explaining 
the origin of high-temperature superconductivity l29l and 
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other emergent states of matter in d- and /-electron systems 
EOl (2H- In addition, spin fluctuations in 5 /-electron systems 
are a manifestation of relativistic effects due to spin-orbit split 
states of order ~500 meV in UCoGas. The emergence of a 
dispersionless band around -80 meV may be related to the 
break-down of the Korringa law for the spin-lattice relaxation 
rate observed above ~75 K (8). In fact, spin-orbit coupling 
may provide a novel electronic tunability for the interaction 
strength and characteristic frequency of the mediating boson. 
The results for the self-energy will motivate efforts to iden- 
tify generic dispersion anomalies and 'waterfall' physics in a 



larger class of actinide and heavy-fermion compounds. 
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VI. SUPPORTING MATERIAL 

VII. MEASUREMENTS 

The photoemission measurements of the electronic struc- 
ture of UCoGas were performed in the Synchrotron Radia- 
tion Center (SRC) using high-quality single crystals grown by 
the flux method. The Plane Grating Monochromator beam- 
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11 momentum points between T and M. Figure Sl(b) shows 
a 3D visualization of the data. The z axis value of the surface 
plot is proportional to the intensity of the spectral function and 
clearly delineates peaks, dips and ridges. The measured data 
points are shown as projected points in the horizontal xy plane 
and as vertices in the overlay ed surface mesh. Figure Sl(c) 
shows a side view of the spectral function to emphasize the 
peak-dip-hump structure seen in the integrated photoemission 
spectrum. The intensity contour map is projected onto the 
horizontal xy plane. 



IX. PHOTON ENERGY DEPENDENCE OF 
HIGH-ENERGY ANOMALY 



FIG. 4. Fig. SI. 3D ARPES visualization, (a) Visualization of the 
experimental ARPES intensity contour map of UCoGas in the Bril- 
louin zone between M and V. (b) Side view of the 3D surface plot 
of the ARPES data with projection of the contour map onto the hori- 
zontal plane, (c) 3D surface plot of the ARPES data with projection 
onto horizontal plane. 



line and a single-channel photoelectron energy analyzer were 
used. The base pressure in the measurement chamber was at 
3 x 10 -11 Torr with cryostat and samples at room temperature 
and at 2x 10 -11 with cryostat and samples at 12 K. The accep- 
tance angle of the analyzer was set to ± 1 degree, and energy 
resolution at 46 eV and 108 eV was 20 meV and 35 meV, re- 
spectively. The EDC were recorded in 12.5 meV increments 
at 1 1 equidistant steps between T and M. The sample was ori- 
ented using a Laue camera prior to the measurement. A de- 
tailed set of normal emission scans was collected in order to 
estimate the inner potential value, determine the correct pho- 
ton energy at the T point (46 eV) in the Brillouin zone, and 
search for k z dispersion. The observed k z dispersion of bands 
rules out the presence of spurious surface states resulting in 
flat bands below the Fermi energy. The difference between 
the 108 eV and 102 eV resonances is due to different cross- 
sections for 5/ and bd electrons, respectively. Hence it is 
characteristic of the strength of the 5/ signal [ 1 ]. The U 5d- 
5f resonance at 108 eV probes the midpoint between T and 
Z with k z ~ 7r/2, while at 102 eV it probes near k z « 0.87T. 
Single crystals of UCoGas and UNiGa 5 were grown from a 
gallium flux l6l . 



VIII. 3D VISUALIZATION OF THE ARPES SPECTRA 

Figure SI shows three different views of the experimen- 
tal ARPES data in the Brillouin zone between T = (0, 0, 0) 
and M= (1,1,0). It facilitates a qualitative analysis of the 
low-energy and high-energy anomalies, kinks, waterfalls and 
dispersionless bands through visual inspection, while at the 
same time providing an overview. The intensity contour map 
is plotted in Fig. SI (a) with small spheres indicating the mea- 
sured points. The EDCs were recorded in steps of 12.5 meV at 



hv = 36 eV hv = 46 eV hv = 56 eV 




FIG. 5. Fig. S2. ARPES spectra at three representative photon 
energies The dashed lines are guide to the eyes about the nature of 
the dispersion. 

Figure S2 demonstrates that the salient high-energy fea- 
ture is a genuine manifestation of the intrinsic correlation ef- 
fect, and is very much independent of the photon-energy. We 
chose these three representative photon energies, yet as wide 
as 20 eV, out of a total of 16 photon energy scans done with 
high resolution ARPES, in which the signature of the high en- 
ergy scale is present in all of them and also they best covers 
one full Brillouin zone (j-T) and midpoint at Z. 

In a separate normal emission study we have sampled over 
30 photon energies between 14 eV and resonances (108eV) 
and the intensity variations in the energy region of interest cor- 
respond to moving through consecutive Brillouin zones, and 
not to the Co 3p to 3d resonances. In 5d-5f resonance study 
we also clearly show that the features at 1 eV binding energy 
are of 5/ character. 

All above findings clearly demonstrate that the high energy 
features we observe are not related to the matrix element ef- 
fects, but are a direct consequence of the electronic structure 
renormalization. We note that the photon energy dependence 
of the 'waterfall' physics obtained earlier for the Bi-based 
cuprate[2] is caused due to the bi-layer splitting! 3], and thus 
not been seen in any other single-layered cuprates,|4] or in 
Sr2Ru204|5 ] or in the present system. 
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A. 



Why is the spectral function anomaly not created by 
ARPES matrix-element effects? 




-1.0 
Energy [eV] 



FIG. 6. Fig. S3. EDC spectra after subtracting the background 
matrix-element contributions. In all figure panels, red solid lines 
are the measured spectra at several representative momentum points. 
The blue dashed lines are the Shirley background for the correspond- 
ing spectra, while green dashed lines are the corrected spectra ob- 
tained after subtracting the background. Even after subtracting the 
background, all spectra continue to exhibit the peak-dip-hump fea- 
ture of our present interest. 



In order to convincingly establish the fact that the spectral 
function anomaly observed in the ARPES spectra is a true 
manifestation of the many-body effects, we adopt a method 
originally proposed by Shirley. Q According to him, the 
ARPES background below the energy scales of the actual va- 
lence bands becomes nearly constant. The background in the 
energy scales of the valence band arises entirely from valence 
band photoelectrons that were inelastically scattered before 
leaving the sample. In this method the 'primary' photoelec- 
trons (un-scattered in the photoemission process) can be sepa- 
rated from the background of inelastic 'secondary' photoelec- 
trons as 1 8 ] 



P{E) = R(E) - B{E) = R(E) -k\ P(E')dE' \ 

Jo 



(3) 



where R(E) is the raw measured spectral function at constant 
momentum or energy distribution curve (EDC), B(E) is the 
background, and P(E) is the primary spectrum or peak of 
the EDC data after Shirley background correction. Here E 
is the binding energy measured relative to the Fermi level at 
E = 0. Eq. (7) is a self-consistent equation, which implies 
that the extrinsic scattering process that a photoelectron un- 
dergoes before being detected depends on its intrinsic spec- 
tral weight distribution. The proportionality constant k deter- 
mines the fraction of primary photoelectrons that participated 
in the scattering process to become the secondary photoelec- 
trons. In principle, k will depend on the band index. However, 



for the energy scale of present interest, we have demonstrated 
in the main text that the bands are well separated up to the 
energy scale of ~ 0.8 eV along M— >• T, and thus we assume 
ft to be independent of energy, but momentum dependent. 

As shown in Fig. S3, the obtained Shirley background term 
is almost featureless at all momenta. The extracted primary 
spectra along the band 'Bl' shown in the main text, as well 
as at three high-symmetry momentum values, continue to ex- 
hibit the peak-dip-hump structure. These results convincingly 
rule out the possibility that such spectral function anomaly is 
extrinsic, but rather a manifestation of intrinsic many-body ef- 
fects. 



X. THEORY 



A. Band structure calculations 




FIG. 7. Fig. S4. GGA band structure. Computed ab-initio band 
structure along high-symmetry line. The results agree well with pre- 
vious calculation using full-potential code. (12) 

We performed electronic structure calculations of UCoGa 5 
within the framework of DFT and the results are shown in 
Fig. S4. Our calculations were carried out by using the 
full-potential linearized augmented plane wave (FP-LAPW) 
method as implemented in the WIEN2k code (9). The gen- 
eralized gradient approximation (GGA) fTOl was used for 
exchange-correlation functional. The spin-orbit coupling was 
included in a second variational way, for which relativistic 
Pi/2 local orbitals were added into the basis set for the de- 
scription of the 6p states of uranium [11]. The energy spread 
to separate the localized valence states was -6 Ryd. The 
muffin-tin radii were: 2.5ao for U, 2.48ao for Co, and 2.2ao 
for Ga, where ao is the Bohr radius. The criterion for the 
number of plane waves was i^^if max = 8 and the number 
of k-points was 40 x 40 x 25. 



B. Spin-fluctuation coupling calculation 

The non-interacting spin susceptibility for a multiband sys- 
tem is a tensor in the orbital basis fT3lL whose components 
are 
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1 

"TV 



k,n,m 



i(fc)^*(fc)C(fc + 9)C(fc + «) 



/(Ufc))~/(em(fc + g)) 

n + ^-(Cn(fe)-Cm(fe + g))" 



(4) 



r 



Here (j) l n is the eigenstate for n th band E n projected on i th - 
orbital, which are obtained directly from first-principles band- 
structure calculations. The susceptibility calculation is done 
with 40 bands that are present in the relevant ±10 eV energy 
window around the Fermi level. The many -body corrections 
are obtained within the random-phase-approximation (RPA) 
given by 



(5) 



The hat symbol A over a quantitity means that it is a matrix of 
N 2 x TV 2 size where N is the number of orbitals. The matrix 
X° has elements Xstpr- The interaction matrix U is defined 
in the same basis consisting of intraorbital U, interorbital V, 
Hund's coupling J and pair-scatterings J' terms (131Q41. In 
the present calculation, we will ignore for simplicity the ma- 
trix elements due to orbital overlap of the eigenstates, i.e., we 
assume (j) l n = 1 when i = n. Such an approximation sim- 
plifies the calculation and x° becomes a diagonal matrix with 
J = J' = V = 0,andZ7=leV. 

In the fluctuation exchange approximation 03), the matrix 
elements of the spin-fluctuation interaction vertex V are given 
by 



Vst pr (q,0) = 



:Ux"(q,W 



(6) 



stpr 



From this expression the self-energy is obtained by using the 
Feynman-Dyson formalism. When employing the spectral 
representation of the Matsubara Green's function matrix G, 



G st (k,iuj m ) = ^ 



E n (k)' 



(7) 



We can write the self-energy within the GW approximation, 
using the fluctuation-exchange interaction matrix V: 



Im£ s t(fc,cj) 



mi: 



procedure is numerically expensive, especially in multiband 
systems. Therefore, we adopt a modified self-consistency 
scheme, where we expand the real part of the self-energy 
E' w (1 — Z~ x )uj in the low-energy region, where E" w 0. 
The resulting self-energy dressed quasiparticle dispersions 
£ n (fc) = Z^ n (k) are used in Eq. (2)-(6), which keeps all the 
formalism unchanged with respect to Z. Finally, it introduces 
a vertex correction which simplifies to qr, co>, Q) = 1/Z 
according to Ward's identity. 



XI. 



DETAILS OF SELF-ENERGY CORRECTION TO 
INDIVIDUAL BANDS 



All bands + E All bands + Re E All bands + Im I 




(f) 














\ '"Ac 











dn [f st (k)+n B (n)}T(k,q, 



pr 

V:}Y(q,n) lmG p 



.(fc — q,uj — Q), 



The above equation denotes the quasiparticle excitation at the 
frequency uj. To reduce the number of indices in the main 
text, we absorb the orbital indices s and t into its correspond- 
ing band index n: T, st —> £ n . The electronic and bosonic 
occupation numbers are f s t(k) and n^(^). The real part 
of the self-energy is obtained from Eq. ^ by employing the 
Kramers-Kronig relationship. Combined with Dyson's equa- 
tion 



G~ l = Gq 1 



(9) 



we obtain the many-body renormalized Green's function G 
matrix. Full self-consistency requires the bare Green's func- 
tion Go to be replaced with self-energy dressed G in x°. This 



FIG. 8. Fig. S5. Computed self-energy corrections to individual 
bands, (a) The self-energy dressed spectral function along M— )► T, 
replotted from Fig. 3(a) of the main text. The Re S [with Im E 
^constant] and Im E [Re E=0] corrected spectral function plots are 
given in (b) and (c). The spectral function renormalization effect due 
to the full self-energy correction is shown individually for various 
ft^nds in (d) through (f). Note that the self-energy itself is calculated 
for all bands and is kept same for all figure panels. 

(8) 

Figure S5 demonstrates how different parts of the self- 
energy dress different electronic states. To make this com- 
parative study easier, we use the same self-energy in all pan- 
els. As shown in Fig. S5(b), the real part of the self-energy 
creates the dispersion renormalization which is strongly en- 
ergy and band dependent. In this case we have used a con- 
stant broadening. The strong band renormalization in the low- 
energy region at M point is created by band 'Bl' [see Fig. S5 
by Re E. By comparing Fig. S5, (b) and (c), we notice that 
although there are renormalized bands present in the energy 
scale of ~ 5 00 meV in 3(b), yet the imaginary part of the self- 
energy creates a dip at this energy as seen in S5(c). This re- 
sult emphasizes the importance of Im E in that it alone can 
produce the peak-dip-hump feature in the spectral function. 
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This is also evident in Fig. S5(e) and S5(f), where we see 
that despite the absence of low-energy quasiparticle bands, 
a nearly dispersionless band is created around -80 meV. For 
a purely momentum-independent self-energy, the low-energy 
feature is completely dispersionless, while momentum depen- 
dence is introduced by the presence of bare electronic bands 
in this energy region, in addition to a momentum-dependent 
self-energy. 

XII. FERMI SURFACE TOPOLOGY RENORMALIZED BY 
INTERACTION AND COMPARISON WITH DE-HAAS VAN 
ALPHEN OSCILLATIONS 



( a ) GGA (b) GGA + 2 dHvA measurements ( c ) 

^^^^^^^^^^^^^^^^ k=n 




r 



FIG. 9. Fig. S6. FS maps and dHvA results, (a) The GGA cal- 
culated FSs for the two bands crossing the Fermi level. The color 
coding reflects the Fermi velocity or the inverse of the density of 
states. The box is centered at the V point, (b) FS cuts at six rep- 
resentative k z -values after including the self-energy (E) correction, 
(c) Measured FS pockets from dHvA oscillations tl6l . In the GGA, 
GGA+E, and dHvA data, all the FS pieces labeled 'a', 'b', and 'c' 
are well reproduced, while the large tube-like FSs 'd' obtained in our 
calculation in (a) and (b) around (7r/2, 7r/2, h) are not seen in the 
dHvA experiment in (c). 

In Fig. S6, we present our theoretical Fermi surface (FS) 
maps calculated without self-energy correction in (a), with 
self-energy correction in (b), and compare with the measured 
FSs via the de-Haas van Alphen (dHvA) measurements by 
Ikeda et al. fT6l in (c). In the k z = plane, we obtain 
three FSs: at the T point (labeled 'a' FS), at (tt, 0, 0) (la- 
beled 'b' FS), and around (7r/2,7r/2, h) connecting the E to 
S points (labeled 'd' FS), and their equivalent points in the 
Brillouin zone. In the k z = tt plane, the 'a' and 'b' FSs disap- 
peared, whereas the 'c' FSs reach maximum size. The GGA 
results are in good agreement with previous electronic struc- 
ture calculations by Opahle et al fl2l . After including the 
spin-fluctuation interaction, the spectral weight across each 
FS varies, but the general shape and position of each FS re- 
mains very much the same. Only the tube-shaped 'd' pockets 
changed from closed to open orbits along the k z direction due 
to the small shift of the chemical potential to lower energies. 
This sensitivity to small changes in the electronic structure 
calculation is consistent with earlier reports of closed versus 
open tube-like pockets fT2l[T7ll . The 'a', 'b' and 'c' pockets, 
seen in the dHvA measurements, are in qualitative agreement 



with electronic structure calculations, as previously discussed 
by Opahle et al. lfT2l . while the 'd' FS was not observed. It is 
possible that the open-loop pocket 'd' was masked by the rela- 
tively broad fundamental frequency of the 'a' orbit. A targeted 
search for the associated dHvA oscillation should resolve the 
so far missing tube-shaped pocket. 

XIII. CONCEPT OF 'WATERFALL' FEATURE AND ITS 
SIMILARITY WITH CUPRATES 




M k T 1)3 l^kp f 



FIG. 10. Fig. S7. Illustration of band renormalization. (a) 

Schematic energy dependence of Re £, containing a sign-reversal at 
a characteristic energy, (b) Im £ is expected to obtain a peak at the 
energy where Re E changes sign, (c) A linear dispersion is shown 
to be renormalized by Re E. (d) Typical form of the spectral weight 
map in the energy and momentum space due to full E -correction, (e) 
ARPES data of UCoGas presented in Fig. 1 of the main text (now 
replotted in different color map), (f) A typical 'waterfall' data of Bi- 
based cuprate superconductor| 18 1 is shown here for comparison. The 
cyan dashed line guides the eyes on how the experimental dispersion 
is deviated from the linear bare dispersion (blue dashed line). 

In Fig. S7, we demonstrate the concept of 'waterfall' 
physics. The basic concept that underlies this lineshape of 
quasiparticles is that the mass renormalization is strongly en- 
ergy dependent and there exists a crossover where the renor- 
malization changes from positive to negative. In this case, the 
strength of the Im £(fc, uo) peak acquires values that split the 
overall quasiparticle states into two energy scales; at low en- 
ergies it induces itinerant quasiparticle states, while at higher 
energies it creates incoherent localized states. These two en- 
ergy scales yield a ubiquitous peak-dip-hump feature, which 
is consistent with photoemission data, see Fig. 1 for a detailed 
comparison. The resulting Re £(fe, uS) acquires an interesting 
energy dependence as shown in Fig. S7(a). The peak height 
of Im E(fe,o;) is sufficiently large such that it leads to Re 
E(fe, uj) at this energy, which is possible within the intermedi- 
ate Coulomb interaction scenario, but not the weak-coupling 
regime. Below this crossover energy Re E(fe, uS) > 0, which 
causes the quasiparticle to acquire a mass enhancement or 
reduction in quasiparticle weight, i.e., Z < 1. As a conse- 
quence the corresponding bands are renormalized toward the 
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Fermi level. Above this characteristic energy Re E(fe, uo) < 0, 
making the quasiparticle states shift toward higher energy 
(Z > 1). At the crossover energy, the dispersion renormal- 
ization is completely eliminated, making the overall quasipar- 
ticle state to sharply drop from Z < 1 to Z > 1 regions. 
At the crossover energy, the spectral weight is strongly sup- 
pressed due to the peak in Im £(fe, uj). The obtained shape of 



the quasiparticle dispersion is dubbed 'waterfall' feature due 
to both its visual and conceptual similarity. 

To demonstrate that the ARPES data of UCoGa 5 indeed ob- 
tain similar 'waterfall' feature as in cuprate superconductors, 
we replot our data in the same colorscale of the data available 
for Bi-based cuprate ifTHl . The characteristic similarity of the 
lineshape of the two spectra affirms our claim that a 'water- 
fall' feature is present in actinide systems. 
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